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Genera]  Description  of  Project 

The  objective  of '‘3 -Dimensional  Crack  Growth  Behavior  of  Turbine  Engine  Materials” 
is  to  evaluate  2-dimensional  and  3 -dimensional  crack  growth  and  material  mechanics  in  turbine 
engine  components  and  its  corresponding  lab  specimens  under  thermomechanical  loading 
configurations.  In  particular,  experiments  are  conducted  on  titanium  and  nickel-basc  super  alloys 
under  simulated  engine  loading  conditions.  The  resulting  residual  stresses,  the  deformation  and 
displacement  of  the  geometry,  and  the  initiation  and  growth  of  cracks  are  analyzed.  Such 
experiments  are  simulated  in  finite  element  analysis  (FEA)  computer  programs.  Previous 
geometries  analyzed  in  this  project  include  turbine  disk  blades,  horseshoe-shaped  beams,  and 
dog  bone  specimens. 

There  are  several  reasons  for  using  FEA  programs.  First,  it  may  be  desired  to  gamer  the 
results  (residual  stress,  crack  growth  behavior)  from  the  FEA  program  and  compare  the  results 
from  the  previous  laboratory  testing  or,  the  given  loading  conditions  may  be  analyzed  in  the 
program  and  the  outcomes  may  be  used  to  predict  the  appropriate  effects  in  the  lab  testing. 
Another  reason  a  program  may  be  used  is  to  determine  the  optimum  geometry  to  avoid  certain 
ranges  of  stresses.  Still  another  need  for  the  analysis  programs  may  be  to  check  the  accuracy  of 
one  FEA  program  to  another;  for  example,  results  may  be  generated  in  both  AN  SYS  (Analysis 
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System)  and  AflAQUS  and  compared  to  theoretical  calculations  to  check  which  program  yields 
better  predictions  for  a  given  problem. 

Description!  of  Research 

Software  was  the  main  tool  in  this  project.  AN  SYS  and  ABAQUS  are  programs  that  can 
model  and  analyse  both  2D  and  3D  problems.  Since  AN  SYS  and  ABAQUS  have  these 
capabilities,  functions  such  as  modeling  and  meshing  a  given  geometry,  adding  appropriate 
material  properties,  creating  analogous  boundary  and  loading  conditions  and  analyzing  the 
problem  are  all  possible  with  both  programs,  For  crack  growth,  FRANC  3D  is  used  to  insert  a 
crack  and  re-mesh  the  present  mesh  in  ANSYS  or  ABAQUS;  FRANC3D  can  also  be  used  to 
extend  existing  cracks  and  create  a  new  mesh.  Then  the  mesh,  with  a  crack  block  included,  is 
analyzed  in  either  ANSYS  or  ABAQUS.  In  conjunction  with  ANSYS  or  ABAQUS,  FRANC3D 
is  also  used  to  compute  stress  intensity  factors. 

A  typical  methodology  is  as  follows,  A  specimen’s  geometry  is  modeled  and  meshed  in 
the  appropriate  program.  Then  material  properties,  such  as  Young’s  Modulus  and  Poisson's 
Ratio,  ate  entered  into  the  program.  The  next  step  is  to  model  the  loading  conditions  under 
which  the  specimen  will  be  exposed  to.  For  example,  typical  loading  of  a  test  specimen  model 
would  be  axial  loading  on  both  ends  with  rigid  body  motion  constraints  mimicking  grips  of  a  test 
frame.  Usually,  the  problem  will  first  be  analyzed  through  these  steps  before  continuing  to  more 
complex  analyses,  such  as  crack  deformation:  this  is  to  ensure  that  the  FEA  program  is 
outputting  appropriate  results  (values  close  to  expected  theoretical  calculations).  Once  the 
problem  is  run  successfully  through  an  FEA  program,  results  are  output. 

The  current  focus  of  analysis  is  the  modeling  and  analysis  of  a  helicopter  lift  frame 
specimen;  see  Figures  1  and  2  for  a  one-quarter  model  of  the  lift  frame.  The  material  is 
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Aluminum  Alloy  7010  T7365  L  The  loads  ate  applied  at  each  bolt  hole  at  each  end  of  the 
specimen,  which  should  result  in  a  nominal  stress  of  about  1  megapascal  (MPa),  Further 
analysis  will  insert  a  crack  along  the  center  hole  of  the  specimen.  The  purpose  is  to  compare 
results  published  in  the  paper  “Life  prediction  from  high  cycle  dynamic  components  using 
damage  tolerance  and  small  threshold  cracks”  by  Vaughan  and  Chang. 

In  AN  SYS,  a  quarter-model  wras  built.  A  one-fourth  geometry  helps  lower  the  analysis 
time  tor  AN  SYS.  Symmetry  conditions  can  be  applied  later  to  the  planes  of  symmetry  to  alert 
the  program  that  the  present  model  is  one-touUh  of  a  total  geometry.  The  following  figures 
showr  two  different  view's  of  the  model. 


Figure  I .  This  is  a  top  view  of  the  specimen  in  ANSYS. 


Figure  2.  This  is  an  isometric  view  of  the  specimen  in  ABAQUS, 

After  the  geometry  was  built,  the  material  properties  for  Aluminum  Alloy  7010  T73651 
had  to  be  entered  in  the  “Materials  Models”  menu  in  ANSYS.  All  model  dimensions  in 
millimeters.  The  Elastic  Modulus  is  entered  as  720,000  N/minA2;  Poisson's  ratio  is  0,33. 
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A  mesh  had  to  be  built  next.  Figure  3  shows  where  hexahedral  and  tetrahedral  element 


types  were  located  in  the  model. 


Figure  3,  The  model  after  a  mesh  has  been  added.  The  boxed  areas  indicate  where  the  elements 
are  tetrahedral;  the  remaining  parts  of  the  mesh  are  hexahedral. 

Next  the  boundary  and  symmetry  conditions  were  added.  The  middle  line  of  the 
specimen  wras  restricted  in  the  X  direction  in  order  to  keep  the  middle  part  of  the  specimen  fixed; 
see  Figures  4  and  5,  Symmetry  conditions  w-ere  added  to  the  bottom  plane  and  to  the  side  plane 
as  indicated  in  Figures  6.  7.  and  8. 


Figure  4.  The  line,  indicated  by  the  three  arrows,  is  where  the  geometry  is  constrained  in  the  x- 
direction. 


Figure  S.  This  is  a  front  view-  to  show  where  the  constrained  line  is  in  relation  to  the  rest  of  the 
model. 
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Figure  6.  The  symmetry  boundary  condition  was  applied  to  the  bottom  and  side  planes,  as 
indicated  here  by  the  boxes. 


Figure  7.  A  view  of  the  bottom  plane  of  symmetry. 


Figure  8,  An  isometric  view  to  show  the  boundary  and  symmetry  conditions. 

To  add  ioads,  two  additional  volumes  were  created  to  mimic  pins  acting  against  the  bolts 
holes.  Figure  9  shows  the  model  with  the  “pins’1  included.  They  were  also  meshed. 
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Figure  9.  The  model  with  the  pins  inserted  into  the  bolt  holes. 

Once  the  pins  were  meshed,  the  necessary  nodal  forces  were  applied  to  each  pin.  Load 
values  such  that  the  resulting  stress  in  the  cross  sectional  areas  of  the  ends  of  the  geometry  were  s 
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1  MPa.  This  is  calculated  to  be  690  N  at  each  bolt  hole.  Since  three  nodes  were  acted  on,  each 


node  underwent  a  force  of  230  N,  shown  in  Figure  1 0. 


Figure  10.  One  of  the  pins  showing  where  the  node  forces  were  applied. 

After  all  of  the  boundary,  symmetry,  and  loading  conditions  were  applied,  the  geometry 
wras  solved  using  AN  SYS. 

Results 

In  the  Post-Processing  menu,  the  x-component  of  stress  wras  the  main  result  of  interest- 
'flu  s  is  because,  along  the  middle  plane,  it  is  the  hoop  stress  for  the  main  hole  in  the  model.  This 
stress  will  cause  the  future  crack  to  open.  Figure  1 1  shows  the  contour  plot. 


,6?1333  E.3S3  4 .  034  5.705 

-.144.333  1 . 5Z7  3,  190  4.:07  6  141 


Figure  1 1 .  The  contour  plot  for  X  stress  after  the  part  has  been  analyzed  in  AN  SYS. 
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Hoop  stress  was  extracted  from  along  the  thickness  of  the  specimen  as  well  as  along  the 
radius;  this  is  where  the  crack  will  be  inserted  later.  See  Figures  12  and  13  to  follow  where  hoop 
stress  wras  extracted. 


Figure  12.  One  path  of  interest  was  along  the  thickness  of  the  main  hole  (the  thick  black  line), 
along  its  line  of  symmetry.  This  path  is  referred  to  as  the  “Mid  Thickness/’ 


Figure  13.  This  is  the  other  path  of  interest  in  the  model.  It  is  along  the  bottom  plane,  along  the 
center  line.  It  is  known  as  the  “Radial  Depth.” 

The  resulting  stress  graphs  are  shown. 

Hoop  Stress  v.  Through  Thickness 


Figure  14.  The  graph  of  hoop  stress  plotted  against  the  through  thickness. 
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Hoop  Stress  v,  Radial  Depth 


Figure  15.  The  graph  of  hoop  stress  plotted  against  the  radial  depth. 

To  check  how  the  ANSYS  results  for  the  highest  hoop  stress  compared  to  theoretical 
values,  the  theoretical  stress  concentration  factor  were  used.  Specifically,  Table  E-15,  Figure  E- 
15-1  from  the  book.  Mechanical  Engineering  Design,  was  referenced. 


0  at  02  0.3  0.4  0.5  O.ti  0.7  0.* 

dtw 

figure  16.  This  is  Figure  E-15-1  from  Table  E-15  in  Mechanical  Engineering  Design , 

The  dimensions  at  the  cross  section  of  the  lift  frame  are  as  follows:  d  (diameter)  =  60mm 
and  w  (width)  -  132mm,  From  these  dimensions,  it  was  found  from  Figure  E-15-1  that  the  stress 
concentration  factor,  K' ,  w^as  about  2,225.  To  find  the  highest  stress  that  occurs  in  the  cross 

section,  K'  is  multiplied  by  the  net  scetionai  stress.  The  net  sectional  stress  is  found  by  the 
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formula  <y,}  =  — ,  where  F  is  the  axial  force  and  A  is  the  cross  sectional  area.  A  is  defined  by  the 


formula  A=  (w-d)  ts  w  being  width,  d  diameter,  and  t  the  equivalent  thickness.  The  following  is 
a  sample  calculation  to  find  the  maximum  theoretical  stress. 


F 


F  =  axial  force 
A  =  (w-d)*t, 

>t  =  width  =  132mm 
d  =  diameter  =  60mm 

x  -  equivalent  thickness  (due  to  the  three  different  thicknesses) 


thickness 

%  width 

sides 

50mm 

0.222222 

11.11111 

majority 

2mm 

0.638369 

1.277778 

hole  lip 

6mm 

0.138889 

0.833333 

From  this  table,  t  =  13.222mm 


Then,  A  =  (132-60)  *  (13. 222)  =  952  mm 2 


F  =  2760N 
F  2760iV 
A  952mm' 


2.899 MPa 


<r0iMt  =K,xcr0=  2225  x  2.899 MPa  =  SASIMPa 


The  net  sectional  stress  was  found  to  be  2,899  MPa;  multiplying  this  stress  by  the  stress 
concentration  factor  gave  us  a  maximum  theoretical  stress  of  6.45 1  MPa.  The  maximum  stress 
found  in  ANSYS  was  6.541  MPa.  This  was  a  1 .407  %  difference  from  the  calculated  theoretical 
value. 


Due  to  the  small  percentage  difference  between  theoretical  and  analytical  maximum 
stress,  it  appears  that  there  is  agreement  between  the  two  results.  The  next  step  in  this  project 
wall  be  to  create  a  crack  block  mesh  and  re-mesh  the  entire  geometry'  using  FRANC3D.  Because 
working  involving  crack  analysis  is  currently  ongoing,  it  is  not  included  in  this  current  technical 
report- 
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